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Four quantities deducible from nuclear structure experiments have been claimed to correlate to 
the slope parameter L of the symmetry energy; the neutron skin thickness, the cross section of 
low-energy dipole (LED) mode, dipole polarizability a_D, and an So (i.e. product of a_D and the 
symmetry energy So)- By the calculations in the Hartree-Fock plus random-phase approximation 
with various effective interactions, we compare the correlations between L and these four quantities. 

The correlation derived from different interactions and the correlation from a class of interactions 
that are identical in the symmetric matter as well as in So are simultaneously examined. These 
two types of correlations may behave differently, as exemplified in the correlation of an to L. It 
is found that the neutron skin thickness and anSo correlate well to L, and therefore are suitable 
for narrowing down the value of L via experiments. The LED emergence and upgrowth makes the 
aoSo-L correlation strong, although these correlations are disarranged when neutron halo appears 
in the ground state. 

PACS numbers: 21.65.-f, 21.65.Mn, 24.30.Cz, 24.30.Gd, 21.60.Jz, 25.20.-x 


I. INTRODUCTION 


Properties of nuclear matter is a basic subject in nu¬ 
clear physics. The equation of state (EoS) of the sym¬ 
metric nuclear matter (SNM), which is characterized by 
the saturation density po, the saturation energy E/A{po) 
and the incompressibility K^o, has been studied for a 
long time and its properties around pq are known rather 
well. In contrast, the EoS of the pure neutron mat¬ 
ter (PNM) has not been established, despite its impor¬ 
tance connected with compact astrophysical objects, e.g. 
neutron stars (NSs). Recent observation of a two-solar- 
mass (2Mq) ns [H has imposed a constraint on the EoS, 
and has given an additional momentum for resolving the 
PNM EoS in particular. Based on the SNM EoS, the 
PNM EoS is mostly governed by the symmetry energy 
S' as a function of density p, which is characterized by 
So = S{p = Po) and the slope parameter. 


L — 3po 


dS{p) 

dp 


( 1 ) 


As So has long been investigated and is known rather 
well, the current uncertainty in the PNM EoS mainly 
originates in the uncertainty in L. 

Although pure neutron many-body systems do not ex¬ 
ist on earth, experiments using radioactive beams dis¬ 
closed that many nuclei have certain volumes dominated 
by neutrons; i.e. neutron skins. This may open a pos¬ 
sibility to constrain the PNM EoS from experiments on 
structure of the neutron-rich nuclei. Objects dominated 
by neutrons may be formed also in the process of nu¬ 
clear reactions, which could leave a signal in observables. 
Many studies narrowing the PNM EoS have been de¬ 
voted to searching observables which strongly correlate 
with L; e. g. nuclear mass systematics , neutron skin 
thickness l7H12| , fragmentation in the heavy ion colli¬ 
sions [l3l - [lq . and low-lying El mode (LED) [13) E] in 


unstable nuclei. Among them, we focus on quantities rel¬ 
evant to structure of specific nuclides, for which model- 
dependence is considered to be relatively weak. 


In Ref. m , the neutron skin thickness Ar. 


np 


208 


Pb 

has been found to correlate linearly to L with a large 
correlation coefficient 0.98, by calculations using 47 ef¬ 
fective interactions. This suggests that accurate deter¬ 
mination of Ar„p serves constraining L. The LED mode 
is considered as a relative oscillation between the neutron 
skin and the remnant core. In Ref. [13 , a linear corre¬ 
lation between the LED cross section (ctled) and L has 
been suggested, from calculations in the random-phase 
approximation (RPA) for ®®Ni and ^^^Sn with 26 effec¬ 
tive interaction. By combining it with the experimen¬ 
tal data, L = 49 — 81MeV has been deduced [T3, 1^ . 
However j the covariance analysis for effective interac¬ 
tions [HHUi has shown that this correlation is not al¬ 
ways strong. Instead, the dipole polarizability ao has 
been claimed to be better in constraining L than cross 
section and transition strength of the LED. If the ao-L 
correlation is assumed, the experimental data in ^°®Pb 
indicate L = 46 ± 15MeV [2^. It has further been ar¬ 
gued, in Ref. [1^, that a product of od and Sq is better 
correlated with L than ud alone, based on the droplet 
model with some assumptions. 


The above four quantities (Ar„p, ctleD) C(d and aoSo) 
have been proposed in separate works, and there have 
been few studies comparing them directly, with excep¬ 
tion of Ref. [1^. Moreover, depending on the studies, 
two different types of the correlations have been argued 
that should be distinguished. The ao-L correlation has 
been investigated using the covariance analysis, for which 
a single interaction and its variants are employed. These 
variants are generated so as to have similar properties 
to the original interaction except L. In contrast, the 
other correlations have been investigated using many in¬ 
teractions with different origin. It is not obvious whether 
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these two types of correlations have the same behav¬ 
ior. We also point out that nucleus-dependence has not 
been discussed sufficiently. Most calculations have been 
implemented in ®®Ni, ^®^Sn and ^°®Pb, partly because 
they are spherical, neutron-rich and accessible by exper¬ 
iments. Nuclear deformation possibly draws complica¬ 
tion, indeed. Still, there could be better candidates. Fur¬ 
ther investigation including careful assessment of correla¬ 
tions is desired in order to constrain L from experimental 
data. 

In this article we investigate the correlations of ^.r^p, 
ctled, o:d and odS'o with L for a number of spherical 
nuclei. The paper is organized as follows. In Sec. m we 
briefly explain interactions we employ and introduce an 
additional term to them, which controls the value of L. 
Numerical results are given in Sec. uni and we discuss the 
interaction- and nucleus-dependence of the correlations. 
Conclusion is given in Sec. IIVI 

II. METHOD 

We perform the RPA calculations on top of the 
Hartree-Fock (HF) wave functions in fully self-consistent 
manner, by using the numerical methods of Refs. [27l.[^. 

In investigating interaction-dependence of the correla¬ 
tions between L and the quantities, we employ a vari¬ 
ety of effective interactions, covering a wide range of L. 
They are three Skyrme interactions which have widely 
been used (SkM* [2^, SLy4 [s^l and SGII ^), two 
latest designed ones (UNEDFO and UNEDFI [32 1), and 
four Skyrme interactions (SkI2, SkI3, SkM and SkI5 [s^) 
that give lar ge L values, and two more Skyrme interac¬ 
tions (SkT4 and Ska M) which are less frequently 
used but useful for checking robustness of the correla¬ 
tions. In addition, three Gogny (D1 [s^, DIS [s^ and 
DIM [^) and two M3Y-type interactions (M3Y-P6 and 
M3Y-P7 [^) are adopted. Using these effective interac¬ 
tions which cover L = 18 — 129 MeV, we discuss the cor¬ 
relations among different interactions (GDI). There have 
been a certain number of relativistic mean-field (RMF) 
calculations. Most of the RMF Lagrangians adopted so 
far tend to give large L values (> 100MeV), which do not 
seem compatible with experimental data. Their results 
are similar, though not identical, to the Skin (n = 2 — 5) 
ones. There may be rooms to obtain RMF Lagrangians 
giving smaller L values. Although we have not imple¬ 
mented the RMF calculations, we shall mention some of 
the RMF results available in literature. 

In the covariance analysis in Refs. (2ll - l^ . a class of 
interactions that share basic properties with an original 
interaction were considered. Following Ref. [i^, we here 
introduce an additional term for the interaction, 

Vij => Vij - Vl [^“(r'i) - Pq ] PaS{r, - Vj ), (2) 

where Pa- is the spin exchange operator. This addi¬ 
tional term does not change So because it vanishes at 
p = po, and has no effects on the SNM EoS because 



FIG. 1: Vl dependence (Eq. @) of (a) EoS and (b) the 
slope parameter L, calculated with SLy4 interaction on set¬ 
ting Vl ~ 0, ±1000, ±2000 MeV, and (c) relation be¬ 

tween neutron skin thickness in ^°®Pb and L shifted by ad¬ 
justing 14- 


{PaS{ri — Vj)) = 0 in the SNM. We thus obtain a 
class of interactions having different L by varying Vl, 
with changing neither SNM EoS nor So- All the non- 
relativistic interactions contain a density-dependent term 
in which the coupling constant is proportional to a power 
of the density. We keep this power a of each original in¬ 
teraction also for the additional term in Eq. ©■ The cor¬ 
relation given by the interactions belonging to the same 
class, which are generated from a single interaction but 
have different Vl, will be called correlation in a single 
class of interactions (GSI) in this paper. 

Figure |T])a) illustrates how Vl affects the EoS, by 
taking the SLy4 interaction and its variants with 
Vl = 0, ±1000, ±2000fm®+®“MeV as an example. 

The L value is changed linearly with Vl, as Vl = 
—2000 (2000) fm®+®“MeV shifts L from the original value 
46 MeV to 17 (75) MeV. The neutron skin thickness is de¬ 
fined by 

= V(r2)„ - ^(r2)p, (3) 

for a specific nuclide. As is expected, L correlates lin¬ 
early with the neutron skin thickness in ^‘^^Pb among 
this class of interactions, as shown in Fig. HJc). The 
additional term changes the binding energy of ^°®Pb by 
~ 15 MeV with Vl = -2000 (2000) fm3+3“MeV. We do 
not take this difference seriously, since this energy shift 
is comparable to the difference of the binding energies 
obtained from different interactions. For instance, UN¬ 
EDFO and UNEDFI yield 1625 and 1643 MeV for ^ospb, 
respectively. 

The El transition operator is expressed as 

q(ei) = ( 4 ) 

i^p iGn 
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after the center of mass correction. Here i is the index 
of nucleons and i G p {i € n) indicates that the sum runs 
over protons (neutrons). The El strength is calculated 


as 




1 


1 


(w - W„) + 7^ (w + OJn) + 7^ 


X (5) 


where n is the index of the excited states and uj denotes 
the excitation energy. For the smearing parameter 7 , we 
adopt 7 = 0.5 MeV, after confirming that the results do 
not change much with 7 = 0.1 — 0.5 MeV. The LED cross 
section ctled is given by 


^^dip 


Ctled — 


9hc 


duJUjS^^^\uj) , 


( 6 ) 


where Wdip is the energy at which is separated 

into the LED and giant dipole resonance (GDR) regions. 
Although the LED and the GDR components could mix 
in certain energy range [dlj, we here separate them by 
energy for simplicity. It is not obvious how Wdip should 
be defined. We determine Wdip as follows. If we find a 
distinguishable LED peak in Wdip is defined as 

the energy corresponding to the minimum of (uj) 
that exists between the LED peak and the GDR. 

The dipole polarizability an is calculated as 


87re2 /■“, 5 (^ 1 ) (a 

an = — 7 — / duj - 

9 7o w 


(7) 
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FIG. 2: (Color online) Correlations of the slope parameter L 
with (a) the neutron skin thickness Ar„p, (b) the LED cross 
section ctled, (c) the dipole polarizability an and (d) an So 
of See text for details. 


III. NUMERICAL RESULTS 
A. Correlations in ^®^Sn 


Owing to the energy denominator, «£> is expected to be 
sensitive to the LED. It should be noted that ao is un¬ 
ambiguously defined unlike ctled- 

As a measure of correlations, it is customary to use the 
correlation coefficient. Eor the two quantities (x, y) for 
which we have data points (xk,yk) (k = 1 , 2 , ,A'd), 

the correlation coefficient is given by 


R[x,y] 


Nd 

'^{xk -x)(yk -y) 



Ni 

Y,ixk-x)\ 

\ k=i \ 

^ivk -y)^ 

fe=i 


( 8 ) 


with x = J2k=i^k/Nd and likewise for y. We obtain 
|i?[a;, y] I = 1 if a: and y are fully correlated and R[x,y] = 0 
if X and y are fully uncorrelated. In the present case 
k corresponds to individual interaction, covering the in¬ 
teractions mentioned above including the variants with 
varying Vl. a; is fixed to be the slope parameter L, and 
y is taken to be Ar„p, ctled, ojd or anSo for a specific 
nuclide. 


EigurelD shows correlations of L with the neutron skin 
thickness Ar„p, the LED cross section ctled, the dipole 
polarizability an and anSo in ^^^Sn, obtained by the 
HF-I-RPA calculations. Effective interactions are distin¬ 
guished by colors and symbols, as listed in the upper part 
of the figure. Results with Vl = 0 are represented by full 
symbols, while those with their Vl ^ 0 variants by open 
symbols. The results of the same class of interactions are 
connected by lines so as to show the GSIs. 

It is seen in Fig. Ufa) that Ar„p correlates well with 
L. Indeed, we obtain R[L, Ar„p] = 0.959. This corre¬ 
lation is well expressed by a linear function, as Ar„p = 
0.00114L -I- O.lBOfm with standard deviation O.OIdfm, 
by assuming the unit of L to be MeV. With respect 
to the CSI, the three interactions UNEDFO, UNEDFl 
and SkI4 give slopes less than 1.0 x 10“^fm/MeV, while 
slopes of the other Skyrme interactions are steeper than 
1.2 X 10“^fm/MeV and those of the Gogny and M3Y 
interactions fall in the narrow range (1.15 ± 0.05) x 
10“^fm/MeV. The maximum (minimum) slope is 1.46 
(0.83) X 10-3fm/MeV of SGII (UNEDFl), which devi¬ 
ates by 30 % from the value htting all the interactions (i.e. 
1.14 X10“^ fm/MeV). We note that slopes of the GSI stay 
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around 1.14 x 10“^ fm/MeV within 10% in more than 
half of the interactions. The GDI (correlations among 
the interactions with Vl = 0) is strong as well, having 
R[L, Ar„p] = 0.939. Thus L can be well constrained 
by Ar„p in ^^^Sn if it is measured precisely. The stan¬ 
dard deviation 0.014 fm is converted to an uncertainty of 
12MeV for L. 

Correlations between L and (Tled are shown in 
Fig. [2Kb). We discard the uled results in the case that 
has two peaks in the LED region, because we 
cannot unambiguously determine Wdip at which the LED 
and GDR regions are separated, and Wdip may change 
discontinuously by changing Vl even if we adopt a cer¬ 
tain definition. Four interactions SkI2, SkI3, SkI4 and 
SkI5 and their variants produce quite large uled, which 
clearly deviate from the results of the other interactions. 
The CSI are not similar even within these four classes of 
interactions. It is noted that the Gogny and M3Y interac¬ 
tions yield correlations similar to the Skyrme interactions 
other than the above Ski series. If we ignore the results of 
the Ski series, i?[L,crLED] = 0.928 is obtained and ctled 
can be fitted to a linear function of L as ctled = 0.399L-I- 
15.4mbMeV with the standard deviation S.OmbMeV. 
When we fit ctled by a quadratic function, we obtain 
ctled = 0.00I38L^-|-0.238L-|-18.7mb MeV with the stan¬ 
dard deviation 4.7 mb MeV. Gompared with Ref. (see 
Fig. 2(b) of Ref. [13)) the slope of the linear function 
is smaller by a factor ^ 2. This discrepancy can be 
interpreted as follows: In Ref. [13) the GDI of ctled 
with L has been investigated via 19 Skyrme interac¬ 
tions and 7 relativistic effective Lagrangians which cover 
L = 0 — 130 MeV. Among them, seven relativistic La¬ 
grangians and three Skyrme interactions SkI2, SkI3 and 
SK255 [13, all of which give L > 100 MeV, seem to be¬ 
have differently from the other interactions. The high 
weight (10 out of the 26 interactions) of these large-L 
interactions leads to the steep slope in Ref. If we 
exclude the results of SkI2, SkI3, SK255 and the RMF 
in Fig. 2(b) of Ref. [13 and refit the others to a linear 
function, the slope is compatible with our result. How¬ 
ever, with ambiguity in the definition of ctled and large 
deviation by certain interactions, we conclude that ctled 
is currently unsuitable for constraining L. 

Figure He) shows the Ud-L relations. Despite the 
relatively large value of R[L,aD\ = 0.90, the lines rep¬ 
resenting the CSI are widely scattered. This indicates 
that the GDI behaves differently from the CSI. If all the 
results of au are fitted to a linear function, we obtain 
ud = 0.026IL-|-5.94e^fm^/MeV with the standard devi¬ 
ation 0.53e^fm^/MeV. However, the slopes given by the 
CSI are significantly larger; 0.031 — 0.051 e^fm^/MeV^ 
with the Skyrme interactions, ~ 0.027 with the Gogny in¬ 
teractions and ^ 0.033 with the M3Y interactions. The 
intercepts are also distributed in as wide range as 2.14 
- 6.15 e^fm^/MeV. It is thus important to take into ac¬ 
count the CSI and the GDI simultaneously. The ao-L 
correlation might look good when we pay attention only 
to the CSI like in the previous covariance analysis, and 



L [MeV] 


FIG. 3: (Color online) Correlations of Ar„p and L in (a) 
^°®Pb, (b) ■‘°’'‘®Ca and (c) ®®’®^Ni. See Fig. E) for colors and 
symbols. 


likewise to the GDI. However, there exists notable differ¬ 
ence between the CSI and the GDI. It is not necessarily 
suitable to constrain L only hy an- 

The aoSo-L correlations are shown in Fig. Hd). The 
strong correlation between an So and L is clearly seen. 
The correlation coefficient is R[L,aDSo] = 0.953. The 
linear fitting gives aoSo = 1.13L -|- 170e^fm^ with 
the standard deviation 15e^fm^, which corresponds to 
13 MeV uncertainty of L, and the quadratic fit gives 
o^dSo = 0.00393L^ -I- 0.617L -I- 180e^fm^ with the stan¬ 
dard deviation 13e^fm^. Even if we restrict ourselves 
to the GDI by setting Vl = 0, the correlations have 
similar behavior; R[L,aL>So] = 0.947, the fitted linear 
function is aoSo = 1.08L -|- 169e^fm^ with the stan¬ 
dard deviation lle^fm^, and the quadratic function is 
aoSo = 0.00208L^ -I- 0.773L -I- 177e^fm^ with the stan¬ 
dard deviation lle^fm^. As pointed out in Ref. [13, So 
has positive correlation with L among the interactions 
with Vl = 0. This helps ajjSo to correlate with L better 
than an alone. Thus aoSo will be useful in constraining 
L, although it requires precise assessment of So- 

We have investigated the correlations of Ar„p, ctleD) 
a LI and ajjSo in ^^^Sn with L, and have found that Ar^p 
and an So are promising for constraining L. 


B. Nucleus-dependence 

The correlations between L and observables related to 
the neutron skin were discussed mainly in ®®Ni, ^^^Sn and 
208pb jjj previous studies. We next consider nucleus- 
dependence of the Avnp-L and aoSo-L correlations. 
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We have calculated the Arnp-L correlations in doubly- 
magic nuclei and in nearly-doubly-magic nuclei, i®’22,24q^ 
40,48,54,70Cg^^ 68.78,84j^-^ 132.140.176 Sj^ 208 pi 3 ^ SOme of 

which are plotted in Fig. [3] In ^°®Pb, Ar„p correlates 
well to L, giving R[L, Ar„p(^°®Pb)] = 0.965. This re¬ 
sult is consistent with that reported in Ref. [HI- The 
linear function obtained by fitting is Ar„p(^°®Pb) = 
0.00107L-l-0.103fm with the standard deviation 0.013 fm, 
being equivalent to 12 MeV uncertainty of L. The slope 
of the fitted function is smaller by ^ 30% than that of 
Ref. 0. This discrepancy is again attributed to con¬ 
tribution of the RMF results with L ^ 100 MeV, because 
they increase the slope in Fig. 3 of Ref. 0. Still the 
Avnp-L correlation in ^°®Pb is so strong to be promis¬ 
ing for getting constraint on L. The Ar„p-L correlation 
gradually becomes the weaker for the lighter nuclei. No¬ 
tice that the steeper slope in the linear function tends 
to make the correlation coefficient the larger. When er¬ 
rors in experimental data are taken into consideration, a 
steep slope is further advantageous in constraining L. As 
mentioned above, we obtain i?[L, Ar„p(^^^Sn)] = 0.959. 
In ®®Ni, i?[L, Ar„p(®®Ni)] = 0.901 and the linear fit¬ 
ting gives Ar„p(68Ni) = 0.000761L-t 0.133fm. In ^SCa, 
the correlation coefficient drops to i?[L, Arjip("^®Ca)] = 
0.785 and the linear fitting results in Ar„p("^®Ca) = 
0.000546L + 0.138 fm. In the Z = N nuclei and 
"‘“Ca, the calculated Ar„p’s are almost independent of 
L. The Avnp-L correlation also becomes weak in drip¬ 
line nuclei such as ®^Ni, as shown in Fig. |31[c). Ar^p is 
strongly affected by the spatial extension of the loosely- 
bound neutron orbits around the neutron Fermi level. In 
nuclei near the neutron drip line, the additional term in¬ 
troduced in Eq. @ with negative Vl, which lowers L, lifts 
up the neutron Fermi level and makes the loosely-bound 
orbits extend significantly. This effect is connected to 
the neutron halo which may irregularly increase Ar„p. 
This mechanism makes the Ar„p-L correlation weaker in 
neutron drip-line nuclei, as is seen in ™Ca and 

i76Sn. 

Therefore, the Ar„p in heavy nuclei distant from the 
drip line may be appropriate in constraining L. Mea¬ 
surement on ^°®Pb seems to provide one of the best pos¬ 
sibilities in this respect. However, despite great efforts 
and progress, it is not yet easy to experimentally deter¬ 
mine Arjip(^°®Pb) with good precision. It should also be 
kept in mind that the Ar„p-L correlation has been inves¬ 
tigated only phenomenologically. Without support from 
quantitatively reliable theories, cross checks from other 
nuclei and/or other quantities are important. 

Let us turn to nucleus-dependence of the anSo-L cor¬ 
relation. Because of the energy denominator in Eq. 0, 
aoSo is rather sensitive to the LED, which emerges and 
grows up beyond the magic numbers N = 14, 28,50 and 
82 1^,1^. We expect that anSo correlates better with L 
as the LED develops in the neutron-rich nuclei. In TableH] 
we list R[L, aoSo] for the stable doubly-magic nuclei and 
neutron-rich nuclei having well-developed LED, 
40.48,54^;^^ 68,84nj^ 208 pi 3 ^ The Optimized 



L [MeV] 


FIG. 4: (Color online) aoSo-L correlations in (a) '*®’®‘^Ca 
and (b) Data of is shifted by 50 e^fm^ to 

accommodate them with the results of ^®^Sn in a single plot. 
See Fig. [2] for colors and symbols. 


values of the coefficients a and b when the calculated re¬ 
sults are fitted as aoSo = aL + b, with the standard 
deviation cr of the fitting, are shown as well. 

The left panel of Fig. |4] illustrates how the LED af¬ 
fects the anSo-L correlation, by comparing the results 
of ^^Ca with those of "‘^Ca. From "‘^Ca to ^'^Ca, a^iS'o 
becomes larger and the slope of the fitted linear func¬ 
tion becomes steeper (0.242 to 0.331 e^fm^/MeV). The 
steep slope is expedient for constraining L from ex¬ 
periment. The LED emergence and development con¬ 
tributes to the aoSo-L correlation. However, in ^'^Ca 
the aoSo-L relation of the M3Y-P6 and P7 interac¬ 
tions deviates significantly from that of the other inter¬ 
actions, while such deviation is not found in ^®Ca. As 
a result, we obtain R[L,aDSo{^'^Ca)] = 0.86, smaller 
than R[L, anSoi'^^Ca.)] — 0.90. This is mainly because 
the M3Y-P6 and P7 interactions produce higher neutron 
Fermi level than the other interactions in ^'^Ca, and gen¬ 
erate the neutron halo when we take Vl < 0. As in 
Avnpi presence of the halo disturbs the correlation, since 
the halo may produce large LED and thereby causes 
large ud- It can be confirmed experimentally whether 
or not ^"‘Ca is a halo nucleus. Suppose that the neu¬ 
tron halo is ruled out, ^'^Ca can be a candidate to con¬ 
strain L from aoSo- Excluding the M3Y interactions, 
we obtain i?[L, Q;£)5'o(®^Ca)] = 0.96 and steeper slope 
(0.40 e^fm^/MeV) in the linear fitting. Also for ®®’®^Ni, 
whereas the slope obtained by the linear fitting becomes 
steeper in ®^Ni, R[L, aD5'o(®"‘Ni)] = 0.66 is small because 
of the neutron halo. Similar trend is seen in the drip-line 
nucleus Before applying auS'o in a certain nucleus 
for constraining L, it should be confirmed that the nu¬ 
cleus does not have halo. 
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TABLE I: Correlation coefficient R[L, anSo], and the optimized values of the coefficients a and b when the calculated results 
are fitted as an So = aL + b, with the standard deviation cr of the fitting. 


nucleus {^Z) 

R[L, aD^o] 

a [e^fmVMeV] 

b 

a 

160 

0.848 

0.068 

9.6 

1.8 

24o 

0.628 

0.058 

9.8 

1.5 

"°Ca 

0.881 

0.210 

33.4 

4.7 

00 

0 

0.898 

0.242 

39.6 

5.0 

®4Ca 

0.857 

0.331 

69.1 

8.4 

®®Ni 

0.929 

0.448 

70.5 

7.6 

00 

0.662 

0.574 

142.7 

27.7 

i“Sn 

0.953 

1.130 

170.1 

15.2 


0.934 

1.354 

213.6 

21.9 

208pb 

0.927 

1.864 

336.8 

31.8 


The correlation coefficients are high both in 
i?[L, Q!DS'o(^^^Sn)] = 0.95 and i?[L, Q!DS'o(^'^°Sn)] = 0.93, 
as presented in the right panel of Fig. 01 Although 
R[L,aDSo] slightly decreases from ^^^Sn to ^^°Sn, the 
slope becomes steeper. Both nuclei are suitable for con¬ 
straining L from an So, if ao is accessible in future ex¬ 
periments. 

The aoSo-L correlation in ^°®Pb has been calculated 
in Ref. employing Skyrme interactions and relativis¬ 
tic Lagrangians, and the linear fitting gives the slope 
a = 2.3e^fm^/MeV and the intercept b = 333e^fm^. 
Compared with our result, the intercept is almost equal 
but the slope is steeper. Another result of the anSo-L 
relation is available from Ref. in which only ao- 
L correlation is calculated with a family of relativistic 
Lagrangian. We can see the aoSo-L correlation us¬ 
ing those results. The fitted linear function represent¬ 
ing the aoSo-L correlation of Ref. has the slope 
a ^ 2.9 e^fm^/MeV and the intercept b ^ 310e^fm^. The 
slope is again steeper than our result while the intercept 
is compatible. Therefore, the currently available RMF 
results increase the slope but have small impact on the 
intercept of the aoSo-L relation. 


C. Comparison with droplet model estimation 

The anSo-L correlation has been suggested in Ref. [1^ 
based on the droplet model under some assumptions. 
The relation of aoSo and L reads 


One of the assumptions in the droplet model is [48|, 0 
9 50 


= ^££4-1/3 ^ 1 
XA- ^qA « 1 , 


( 10 ) 


where Q is the surface stiffness coefficient connected with 
the nuclear surface symmetry energy 0 . For the droplet 
model estimation (Eq. ©) to be justified, xa should be 
sufficiently small. In evaluating xa, we use an approxi¬ 
mate expression for Q 



( 11 ) 


instead of calculating Q^in the asymmetric semi-infinite 
nuclear matter P. l48l - [50| . Here a (~ 0.55 fm) is the dif¬ 
fuseness of the symmetric semi-infinite nuclear matter 0 
and Ksym is the 2nd derivative of the symmetry energy 
with respect to the density at the saturation point. The 
calculated xa values are listed in Table HIl accompanying 
Pq, K OD : So, L, RTgyjj^ and Q. 

Figure [S] shows the ratio of calculated with 

HF-bRPA to (ai55o)DM in '^^Ca, ®®Ni, and ^ospb. 
The ratios clearly deviate from unity even for small 
Xa ■ Although we have phenomenologically confirmed the 
anSo-L correlation, the droplet model is not necessarily 
appropriate for justifying the correlation. 


IV. CONCLUSION 


7re 


(au^ojDM ^ ^ 


1 + 


5 L po- PA 


9 So po 


(9) 


where pA 0.1fm“^ [3, While we have phenomeno¬ 
logically confirmed the anSo-L correlation in preceding 
sections, it deserves investigating validity of this rela¬ 
tion. For a given interaction and a nucleus we evaluate 
(a£)5o)DM from Eq. ©, and compare them to aijSo ob¬ 
tained from the HF-I-RPA calculations. 


We have investigated the correlations of L with the fol¬ 
lowing four quantities; the neutron skin thickness Ar„p, 
the cross section of the low-energy dipole (LED) mode 
ctled, the dipole polarizability ao, and the product of 
ao and the symmetry energy So- In order to directly 
compare them and to unravel disorder in observables con¬ 
straining L, we have simultaneously discussed the corre¬ 
lations derived from different interactions (GDI) and the 
correlation in a single class of interactions (CSI). For the 
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TABLE II: Saturation density po, incompressibility of symmetry nuclear matter Kao, symmetry energy So, slope parameter L, 
incompressibility of symmetry energy TsTsym, surface stiffness parameter Q and xa multiplied by the mass dependence, given 
by the Skyrme, Gogny and M3Y interactions with Vl = 0. 



Po [fm“®] 

Kao [MeV] 

So [MeV] 

L [MeV] 

K„y^ [MeV] 

Q [MeV] 

Xa X AU® 

SkM* 

0.160 

216.4 

30.0 

45.8 

-155.8 

35.9 

1.88 

SLy4 

0.160 

229.9 

32.0 

45.9 

-119.7 

42.8 

1.68 

SGII 

0.158 

214.5 

26.8 

37.7 

-145.8 

33.9 

1.78 

UNEDFO 

0.160 

229.8 

30.5 

45.1 

-189.6 

35.8 

1.92 

UNEDFl 

0.159 

219.8 

29.0 

40.0 

-179.4 

35.8 

1.82 

SkI2 

0.157 

240.7 

33.4 

104.3 

70.6 

26.6 

2.82 

SkI3 

0.158 

258.0 

34.8 

100.5 

72.9 

30.2 

2.60 

SkI4 

0.160 

247.7 

29.5 

60.4 

-40.6 

31.9 

2.08 

SkI5 

0.156 

255.6 

36.6 

129.3 

159.4 

27.3 

3.02 

SkT4 

0.159 

262.9 

35.5 

94.1 

-24.5 

30.6 

2.60 

Ska 

0.155 

235.3 

32.9 

74.6 

-78.4 

31.5 

2.35 

D1 

0.166 

229.4 

30.7 

18.4 

-274.6 

41.4 

1.67 

DIS 

0.163 

202.9 

31.1 

22.4 

-241.5 

41.0 

1.71 

DIM 

0.165 

225.0 

28.6 

24.8 

-133.2 

41.0 

1.57 

M3Y-P6 

0.163 

239.7 

32.1 

44.6 

-165.3 

31.8 

2.27 

M3Y-P7 

0.163 

254.7 

31.7 

51.5 

-127.8 

29.2 

2.45 


1.4 H^-'-'---'- 



0.4 H-1-1-1-1-h -j-1-1-1-1- 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Aa 


FIG. 5: (Golor online) Ratio of an So obtained from the 
HF+RPA calculations to the droplet model estimate in ®‘*Ga, 
®®Ni, ^®^Sn and ^®®Pb. See Fig. for colors and symbols. 


latter we introduce an additional term to each interac¬ 
tion, which enables us to control the value of L without 
influencing SNM EoS and Sq. 

The Avnp correlates almost linearly with L in 
heavy nuclei, although there remains slight interaction- 
dependence as recognized via comparison with the re¬ 
sults in Ref. [ll|. The ctled-Tj correlation has a signif¬ 


icant interaction-dependence. Together with ambiguity 
in its definition, uled is not recommended to constrain¬ 
ing L. In the ao-L correlation, we have found that the 
CSI and the GDI behave differently. It is not reasonable 
to constrain L only from an- The aoSg-L correlation 
works well for narrowing down L. The Ar„p and aoSo 
are promising for constraining L, though with ~ 12MeV 
uncertainty. 

The nucleus-dependence of the Avnp-L and aoSo-L 
correlations has also been discussed. While the neu¬ 
tron halo makes the correlations weak, these correlations 
are strong in neutron-rich medium- or heavy-mass nu¬ 
clei without neutron halo. Except neutron-halo nuclei, 
the LED makes the aoSo-L correlation strong and the 
slope of the linear function steep, to which the HF-I-RPA 
results are well fitted. Consequently, the neutron-rich 
nuclei having well-developed LED (e.g. ^"*’Ca and ^'^°Sn) 
are good candidates for obtaining the constraint on L, as 
well as the doubly magic nuclei ^^^Sn and ^°®Pb. 
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